Non-equilibrium two-phase coexistence in a confined granular layer 
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We report the observation of the homogenous nucleation of crystals in a dense layer of steel spheres 
confined between two horizontal plates vibrated vertically. Above a critical vibration amplitude, 
two-layer crystals with square symmetry were found to coexist in steady state with a surrounding 
granular liquid. By analogy to equilibrium hard sphere systems, the phase behavior may be explained 
through entropy maximization. However, dramatic non-equilibrium effects are present, including a 
significant difference in the granular temperatures of the two phases. 

PACS numbers: 45.70.-n,05.70.Fh,05.70.Ln,83.10.Rs 



Statistical mechanics provides a powerful formalism for 
predicting the behavior of systems at or near equilibrium. 
Many natural phenomena, however, occur far from equi- 
librium, and extensions of the machinery of statistical 
mechanics to situations where significant energy flows are 
present would have a wide range of potential applications. 
Some success has been achieved in extending statistical 
mechanics to far-from-equilibrium situations. For exam- 
ple, techniques for calculating the relative probabilities 
of different configurations have been successfully devel- 
oped for non-equilibrium steady states in a few restricted 
situations [jj, I4 yj- Also, effective temperatures based 
on the fluctuations in non-equilibrium steady states have 
been developed recently and, in some cases, have been 
found to equilibrate across different fluctuating quanti- 
ties HUlHSHHEil- To generalize these initial 
successes to a broader spectrum of phenomena and to 
develop insights into other facets of a statistical theory 
of far-from-equilibrium systems, new model systems for 
investigating non-equilibrium steady states must be de- 
veloped and studied. Here we report our investigations 
of a simple, far-from-equilibrium granular system that 
shows that some mechanisms that operate in equilibrium 
appear to persist into situations far-from-equilibrium, 
whereas other basic tenets of equilibrium statistical me- 
chanics must be substantially modified. 

Granular materials are ubiquitous in nature and show 
a remarkable range of non-equilibrium behavior [lit H^ . 
ll3T |. Dynamic steady-states, achieved when energy input 
from an external source balances energy lost through in- 
elastic collisions, provide an ideal testing ground for ex- 
tensions of equilibrium statistical mechanics 0, 0, 0, El ■ 
In this Letter, we report experimental measurements and 
computer simulations of the dynamics of spherical par- 
ticles confined between two horizontal vibrating paral- 
lel plates. We observe a freezing transition from a ho- 
mogeneous, disordered liquid to an ordered solid with 
square symmetry coexisting with a surrounding liquid. 
An essentially identical transition is observed in con- 
fined hard-sphere colloidal suspensions in equilibrium 
jl4l llfil llfil Il7l Il8| , where it is driven by entropy max- 
imization. The presence in both the granular and the 



colloidal system of a solid phase with the same unex- 
pected symmetry which occurs under the same geomet- 
ric constraints and at similar densities strongly suggests 
a common mechanism. Unlike the equilibrium system, 
however, we find that the coexisting phases have dramati- 
cally different granular temperatures, demonstrating that 
the "zeroth law" of thermodynamics is not followed by 
the granular temperature. Taken together, these results 
show that the driving mechanism behind an equilibrium 
phase transition may still operate far from equilibrium, 
but that a thermodynamic theory must account for the 
absence of equipartition in the kinetic energy of the par- 
ticles. 

The granular system is sketched in Fig. Ufa). Pre- 
vious studies which used a similar geometry at lower 
shaking amplitude and lower density than what is pre- 
sented here found a range of complex non-equilibrium 
phenomena including inelastic collapse |l9l l2fj| . hexago- 
nal ordering j 1 OL 12 1| - non-Gaussian velocity distributions 
jlfll Efl E3 | , and velocity correlations [2~j| ■ In our exper- 
iment, stainless steel spheres of diameter a = 1.59 mm 
were placed between a smooth anodized aluminum plate 
and an 11 mm thick Plexiglas lid with an anti-static coat- 
ing. A gap spacing of 1.75 a between the plate and the 
lid was maintained by circular rings of aluminum and 
Mylar spacers. Using an electromagnetic shaker, the sys- 
tem was driven sinusoidally in the vertical direction with 
frequency v and amplitude A. The motion of the balls 
was imaged from above using a high resolution camera 
(Pulnix TM1040). The results presented here were ob- 
tained for densities p = N/N max ranging from 0.8 to 
0.9, where N is the number of balls in the system, and 
Nmax = 11377 is the maximum number of balls that can 
fit in a single hexagonally close-packed layer at rest on 
the bottom plate. For modest vibration amplitude, the 
system appears liquid-like. As the vibration is increased, 
small independent unstable crystalline structures form. 
Increasing the vibration amplitude further causes the 
crystallites to become stable and eventually coalesce to 
form a single crystal which coexists with a surrounding 
granular liquid (see Fig. QJb) and Supplemental Movie 
1 [2i| ) The crystal consists of two layers, each with a 
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FIG. 1: Two phase coexistence in steady state, (a) Side view 
of experimental setup. The system is shaken vertically and 
imaged from above, (b) Experiment: time-averaged image 
of ball positions. Only the top layer of the crystal is visible 
and there are several vacancies (A — 0.085 a, p — 0.9, v — 
80 Hz, averaged over 1 second.) See also Supplemental Movie 
1 |24|. (c, d) Simulation: 3D rendering of instantaneous ball 
positions. In (c), balls in the crystal are colored red, balls in 
the liquid are colored blue. In (d) a close-up of the crystal is 
shown with the top layer transparent. (A — 0.13 a, p = 0.89, 
v = 60 Hz, and 5,000 balls.) 



square symmetry. The balls in the second layer are above 
the centers of the squares formed by the balls in the bot- 
tom layer. The crystals are not densely packed and the 
balls constantly jitter around in the cage formed by their 
neighbors and the confining plates. Rearrangements con- 
tinually occur at the interface, but the average size of the 
crystal does not change as long as the shaking amplitude 
and frequency are held constant. 

In order to verify that the coexistence was not due to 
non-uniformities in the experimental apparatus and to 



measure quantities not readily accessible in the experi- 
ment, we performed molecular dynamics simulations us- 
ing a model that has accurately reproduced many of the 
phenomena observed in a similar system [23l |25| . Ball- 
ball, ball-plate, and ball-lid interactions are characterized 
by three forces: an elastic restoring force, a dissipative 
normal force which produces a velocity-independent co- 
efficient of restitution, and a dissipative tangential fric- 
tion. Periodic boundary conditions in the horizontal 
plane were used. The simulations reproduced the phase 
coexistence (Fig. [flc, °0) an d au of the general features 
of the observed phenomena, such as the existence of a 
critical threshold to nucleation and evaporation. 

Perhaps the most surprising aspect of the transition 
we observe, the presence of a square symmetry instead of 
the hexagonal ordering naively expected for hard-sphere 
interactions, appears to be closely related to the phase 
behavior of solutions of hard sphere colloidal particles 
at simil arly high densities in similar confining geometries 
|l4 Il5l llo ] . For hard spheres the equilibrium configura- 
tion is determined by entropy maximization. For a range 
of gap spacings, including the spacing used in our system, 
two square layers pack more efficiently than two hexago- 
nal layers, thereby maximizing the free volume available 
for each particle and therefore the entropy of the sys- 
tem 0, . The observation of a transition that closely 
matches an equilibrium, entropy driven phase transition 
suggests that a generalized free energy functional might 
be found which describes the behavior of some driven 
granular materials. 

Equilibrium two-phase coexistence requires that the 
two phases have equal temperatures, pressures (apart 
from surface tension corrections) and chemical poten- 
tials. Recent work has focused on extending the concept 
of the chemical potential to non-equilibrium coexistence 
in which the first two conditions are satisfied [2(| W\ . 
To test whether these two conditions are met in this sys- 
tem, we measured the granular temperature T g =< vf >, 
where Vj is a horizontal component of the rapidly fluctu- 
ating velocity of a particle. In the experiment, particle 
displacements were measured using the PIV technique 
described in Ref. [2j|. Using the method described 
in Ref. (2^, we verified that the temporal resolution 
was sufficient to accurately measure instantaneous veloc- 
ities. In both the experiments and simulations, we inves- 
tigated whether the granular temperature equilibrated 
to the same value in the two phases. As shown in Figs. 
Ufa) and [2Jb) , T g is dramatically lower in the crystal 
than in the surrounding liquid, both in the experiments 
and in the simulations. The spontaneous separation into 
phases of different temperatures in a homogeneous sys- 
tem of identical particles is a striking effect that will 
have to be incorporated into models of non-equilibrium 
phase coexistence. It is somewhat reminiscent of 'inelas- 
tic collapse' 19], but in that case the absence of any sig- 
nificant granular temperature in the solid phase arises 




FIG. 2: Temperature and pressure fields, (a) Temperature 
field near the liquid-solid interface, measured in the experi- 
ment, averaged over about 40 seconds. The crystal is in the 
upper left, (p = 0.85, v — 60 Hz, A = 0.154 a, displacements 
measured over 0.33 ms.) (b), Temperature field and (c) pres- 
sure field from the simulation shown in Fig. lc, averaged 
over 2.5 seconds. The crystal is in the center. The pressure 
was calculated using the virial expression, as described in Ref. 
M. 



from the bistability of the ball-plate dynamics at low vi- 
bration amplitudes [2^, |2l| . The results described here 
are observed at high vibration amplitudes where there is 
continuous energy input from the plate into both coexist- 
ing phases. The pressures of the two phases calculated in 
the simulations have nearly the same value, but is slightly 
smaller in the solid phase (Fig. GJc)). 

To further study the properties of the phase coexis- 
tence, we investigated the nucleation of the crystalline 
phase. Starting with the vibration amplitude at a low 
value, we slowly increased the intensity of shaking and 
measured the amplitude at which the crystalline phase 
first nucleates. This procedure was repeated for several 
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FIG. 3: (a) Frequency dependence of the critical amplitudes 
Ac for nucleation (open circles) and evaporation (black disks) 
at p = 0.85. (b) High frequency plateau value of Ac as a 
function of p. (c) Average number of balls of the crystal, No, 
as a function of p (u = 60 Hz, A = 0.145 a). Each point 
represents the average of 10 measurements well separated in 
time . The solid line is a linear fit to the points. 



densities between p = 0.8 and p = 0.9, and for frequen- 
cies between v — 45 Hz and v = 100 Hz. A typical curve 
of the frequency dependence is displayed in Fig. |3[a) for 
p = 0.85. In addition to this 'nucleation line', we also de- 
termined an 'evaporation line' by slowly decreasing the 
amplitude until the crystal disappeared. For v greater 
than 60 Hz, we found that the critical amplitudes were 
roughly independent of frequency. This high frequency 
behavior was found at all densities, but the cutoff frc- 
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quency increased as the density increased. One possible 
explanation for a frequency-independent critical ampli- 
tude is that the vibration may effectively compresses the 
layer. If the balls are moving slowly compared to the 
plate and lid, then they will be mostly confined between 
the maximum plate height and minimum lid height. This 
increase in the density of the system favors nucleation of 
the crystal. This frequency-independent behavior can- 
not persist to low frequencies, however, because the ac- 
celeration, which is proportional to v 2 , must be signifi- 
cantly larger than that due to gravity for the balls to have 
enough kinetic energy to reach the second layer. We used 
the average value of the amplitude in the high frequency 
plateau (60-100 Hz) to define a critical amplitude, Ac, 
and we constructed a "phase diagram" of Ac versus p 
(Fig. EJb)). Ac varies from roughly 0.06 a at p — 0.9 to 
0.1 o at p = 0.8. 

We measured the dependence of the crystal size on the 
number of particles in the system by analyzing the im- 
ages to extract Nc, the number of spheres in the crystal 
|31| . As shown in Fig. Ofc), Nc varies linearly with p 
and extrapolates to zero at p = 0.78. In steady state co- 
existence, the edge of the crystal is "in equilibrium" with 
the surrounding liquid of density p^. Assuming that the 
densities of the coexisting crystal and liquid are inde- 
pendent of the size of the crystal, the area A occupied 
by the crystal of density pc should satisfy the relation 
A/ At = (p — Pl)/{pc — Pl) where At is the total sur- 
face area of the plate, so that Nc oc (p — pl)- The value 
of pl found by extrapolating to Nc = agrees with di- 
rect measurements of the density of the granular liquid 
in the coexistence region. 

No formalism exists for incorporating the entropy into 
a predictive theory on non-equilibrium phase transitions, 
but our results indicate which parts of the equilibrium 
framework need modification. The large difference in the 
granular temperature of the coexisting phases demon- 
strates that the "zeroth law" of thermodynamics is not 
satisfied by the granular temperature. An effective tem- 
perature that does meet this requirement is probably a 
necessary ingredient of a quantitative theory of the phase 
coexistence. By comparing the system described in this 
Letter with the analogous and well understood equilib- 
rium system, new approaches for incorporating the effects 
of forcing and dissipation into a statistical mechanics of 
non-equilibrium phase transitions can be developed and 
tested. 
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